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HONGYU LIU AND TING ZHOU 

Abstract. We investigate two-dimensional invisibility cloaking via trans- 
formation optics approach. The cloaking media possess much more sin- 
gular parameters than those having been considered for three-dimensional 
cloaking in literature. Finite energy solutions for these cloaking devices 
are studied in appropriate weighted Sobolev spaces. We derive some 
crucial properties of the singularly weighted Sobolev spaces. The invis- 
ibility cloaking is then justified by decoupling the underlying singular 
PDEs into one problem in the cloaked region and the other one in the 
cloaking layer. We derive some completely novel characterizations of the 
finite energy solutions corresponding to the singular cloaking problems. 
Particularly, some 'hidden' boundary conditions on the cloaking inter- 
face are shown for the first time. We present our study for a very gen- 
eral system of PDEs, where the Helmholtz equation underlying acoustic 
cloaking is included as a special case. 



1. Introduction 

A region is said to be cloaked if its contents together with the cloak 
are indistinguishable from background space to certain exterior detections. 
Blueprints for making objects invisible to electromagnetic waves were pro- 
posed by Pendry et al. [TT] and Leonhardt [11] in 2006. In the case of 
electrostatics, the same idea was discussed by Greenleaf et al. [5] in 2003. 
The works [5J [11] [T7] rely on transformation optics for the construction of 
cloaking devices, which we shall further examine in the present paper. For 
state-of-the-art surveys on the rapidly growing literature and many applica- 
tions of transformation optics, we refer to [H 1161 H5]- 

In transformation optics, the key ingredient is that optical material pa- 
rameters have transformation properties which could be pushed forward to 
form new material parameters. Then, to construct cloaking devices, the idea 
is to blow up a point in the background space to form the cloaked region. 
The ambient background medium is then pushed forward to form the cloak- 
ing medium. Since the blowing-up transformation is singular, the resulting 
cloaking medium is inevitably singular. Two theoretical approaches have 
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come out to handle the singular cloaking problems. In Kohn et al. [10j . 
the authors introduce the notion of near-invisibility cloaking in electrostat- 
ics from a regularization viewpoint. The singular ideal cloaking is regarded 
as the limit of the regular near-invisibility cloaking depending on certain 
regularizer. For acoustic cloaking, the near-invisibility is investigated in 
j9j [T2l [15] in both two and three dimensions. In [3], Greenleaf et al. pro- 
posed to investigate the physically meaningful solutions, i.e. finite energy 
solutions, corresponding to degenerate differential equations underlying the 
three-dimensional cloaking. The proposal has been shown to work for both 
acoustic and electromagnetic cloaking, and can treat cloaking of passive ob- 
jects as well as active/radiating objects. On the other hand, the analysis 
in [3] is conducted in the geometric setting by taking advantage of the one- 
to-one correspondence in R 3 between the optical parameters and a smooth 
Riemannian metric. This argument does not carry over to M 2 . As one shall 
see in Section 2, for two-dimensional cloaking, the cloaking medium has both 
degeneracy and blow-up singularities at the cloaking interface, making the 
problem more difficult to analyze. 

In this paper, we consider the two-dimensional invisibility cloaking for 
a very general system of second order partial differential equations. The 
Helmholtz equation underlying the acoustic cloaking is included as a spe- 
cial case. In order to handle the singular cloaking problem, we follow the 
finite energy solutions approach from [3]. Recently, Hetmaniuk and Liu 
[6] introduce weighted Sobolev spaces with degenerate weights for three- 
dimensional acoustic cloaking problems, encompassing and generalizing the 
idea of finite energy solutions approach. For the present two-dimensional 
cloaking problems, we study weighted Sobolev spaces with more severely 
singular weights. The invisibility justification follows by study of weak solu- 
tions from the introduced weighted Sobolev space to the underlying singular 
PDEs. Our analysis is given in a very general setting. The cloaking is shown 
to work for arbitrary positive frequency and can cloak both passive media 
and source/sink inside the cloaked region. Since the cloaking media are 
more singular than those for the three-dimensional cloaking, we derive com- 
pletely different and novel characterizations of solutions to the underlying 
wave equations compared to those derived for the three-dimensional acous- 
tic cloaking in [31[6]. Moreover, some 'hidden' boundary conditions on the 
cloaking interface are shown for the first time, giving more insights into the 
invisibility cloaking. 

In this paper, we focus entirely on transformation-optics-approach in con- 
structing cloaking devices. But we mention in passing the other promising 
cloaking schemes including the one based on anomalous localized resonance 
[13] . and another one based on special (object-dependent) coatings [JJ. The 
rest of the paper is organized as follows. In Section 2, we introduce the 
transformation optics and invisibility cloaking and give the construction 
of the two-dimensional radial cloaking devices. Section 3 is devoted to the 
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analysis of the radial cloaking by considering finite energy solutions in singu- 
larly weighted Sobolev spaces. In Section 4, we extend our study to general 
invisibility cloaking. 



2. TRANSFORMATION OPTICS AND RADIAL INVISIBILITY CLOAKING 

We first fix notations for some function spaces which are crucial for our 
study. Let O be a bounded Lipshitz domain in M 2 . Let m > 1 be an 
integer, and u : £1 — > C m be a complex vector-valued function. L p (£l) m is 
the space consisting of C m -valued measurable functions whose components 
belong to D>(Q). Following Schwartz, put g (n) m = C°°(12) m and 3>{n) m = 
C^ mp (Q) m be complex-valued, smooth test function spaces. Let H S (Q) 
denote the standard Sobolev space of order s and 

H S (Q) = closure of in H s {W l ). 

Note that H~ S (Q) is an isometric realization of H S (Q,)* for s G E. The 
definition of the vector Sobolev spaces on 0, H s (n) m = H s (£l;C m ) and 
H s (n) m = H s (tl;C m ) etc. shall now be obvious. 

Next, we introduce the second-order partial differential operator (PDO) 
V of the form 

(2.1) ^-j^^-^o^, 
where the coefficients 

(2.2) A al3 = [0$] G L°°(n) mxm , B = [bpq] G L°°(^) mxm 

are functions from £1 into C mxm , the space of complex mxm matrices. Here 
u) G R and m > 1 is an integer and thus, V acts on a (column) vector-valued 
function u : Q — > C m to give a vector- valued function Vu : ft — > C m , whose 
components are 

m 2 m 

(vu) p = - Yl da (%k d P u k) -^Yl b Pi u i- 

fe=la,/3=l q=l 

In the sequel, let A := [A a/3 ] G C 2mx2m be the block matrix, and V[a,b] be 
the PDO (|2.ip associated with A and B. For the present study, we always 
assume that 

(2.3) {A a P)* = A? a for a, = 1, 2, 

where * denotes the conjugate transpose of a matrix or vector. Moreover, 
we introduce the following algebraic conditions for the coefficient matrices 
A a @ and B. Let < c\ < < oo be two constants. For all x G fi, and 
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arbitrary ^1,^2 G C m and 7/ G C m , we have 

22 2 

(2.4) Cl I6i 2 < » E i^V^m < c 2 x; i&i 2 , 

i=l a,/9=l «=1 

(2.5) cilr?! 2 < ^[B7]]*ri < c 2 \r]\ 2 . 

Next, we associate V\a,b] with a sesquilinear form Q[a,b]i defined by 

(2-6) Q [A ,B](u,v) = [ ( V (A^^«)*a a i; - ) dx. 

Now, we are ready to present the PDE system for our study, 

(2.7) V\a,b]U = f on n, 

where / G if~ 1 (J7) m . The Dirichlet-to-Neumann (DtN) map 

(2.8) a a,bj- H l ' 2 {dtt) m -^H^ 2 {dny n , 

associated with (|2.7|) is defined by 

(2.9) A A,Bj(h) = Yl ^A a ^{d p u) on dfl, 

o,/3<2 

where v = {^0)^=1 is the outward unit normal of d^l, u G H 1 ^) 771 solves 
(|2.7p with = h and 7 is the trace operator for fi. The weak solution in 
(12.91) is variationally given by 



(2.10) Q[a,b](u,<p) = (f,<p)n ■= / /V dx. 



Due to (|2.4j) and (|2.5p . we know that Q[a,b] (and so is "P^s]) is coercive on 
i7 1 (fi) m in the sense that 

(2.11) m(u,u) > ci||u||^i (n)m - c 2 \\uf L 2 (n) m for u G H 1 ^)™. 

Hence, (|2.10p is uniquely solvable except at a discrete set of eigenvalues for 
uj 2 (see, e.g. [13]). In fact, we know that B j is a well-defined continuous 
and invertible operator provided (I2.4|) and ()2.5p are satisfied and u 2 avoids 
the (discrete set of) eigenvalues (cf. [13]). 

In the case m = 1 and yl, £> are both real, (I2.7P is the scalar Helmholtz 
equation. It describes the time-harmonic solutions p{x) = u(x)e~ tut of the 
scalar wave equation pu — B^ 1 \7 ■ (AX7p) = f(x)e~ iljt . Here / represents a 
source/sink inside the region f2. A and B are the acoustic material param- 
eters of the medium supported in Q, related respectively to, density tensor 
and modulus. For a regular acoustic medium, (12. 4p and (|2.5p are the physical 
condition on the material parameters. According to our earlier discussion, 
these are also mathematical conditions to guarantee the well-posedness of 
the underlying Helmholtz equation. In the sequel, we let {£l;A,B,f} de- 
note the medium and the source/sink supported in f2. This is the prototype 
problem of our present study. Moreover, we are concerned with the inverse 
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problems of identifying the inside object {f2; A, B, /} by the exterior wave 
measurements, which are encoded into the DtN operator (|2.8p . The inverse 
problems have widespread practical applications in science and engineering, 
and have received extensive and intensive investigations in last years (see, 
e.g., [EJUS]). In this context, an invisibility cloaking device is introduced as 
follows (see also [9] ) . 

Definition 2.1. Let Q and D be bounded domains in ]R 2 with D ^ Q. Q\D 
and D represent, respectively, the cloaking region and the cloaked region. 
{Q\D; A c , B c , f c } is said to be an invisibility cloaking device for the region 
D with respect to the regular reference/background space {Q; Ab, Bb, /{,} if 

*A M = ^A M fOT a11 " > °> 

where the extended medium {Q; A e , B e } and the extended source f e are 
given by 

({n\D;A c ,B c ,f c } in Q\D, 
\{D;A a ,B a ,f a } in D, 

with {D; A a , B a , f a } arbitrary but regular. 



{O; A e , B e , / e } 



In Definition 12.11 {£1;A,B} is regular means that A and B satisfy the al- 
gebraic conditions (|2.4p and (I2.5p . According to Definition 12. 1\ we note 
that the cloaking device {ft\D; A c , B c , f c } makes the interior target ob- 
ject {D; A a , B a , f a } indistinguishable from the reference/background space 
{f2; Ab, Bf>, f } by exterior detections. In fact, the exterior observer would 
not even be aware that something is being hidden. 

Next, we present the transformation properties of material parameters, 
which are the cruxes of the construction of cloaking devices via transforma- 
tion optics approach. In the following, we let y = F(x) be a diffeomorphism 
such that F : £1 —> Q. Then, the push-forwards of material parameters are 
given by 

F*{fi; A, B} = {n;F*A,F*B} := {n ; A,B}, 



with 



dy P dy„ 



, p,q = 1,2, 



(2.12) A«(y) = V ^L^LA«P{x)J^ 

(2.13) B(y) = B{x).r\ =F -, {y) , 

where J := det([dy a / dxp]) , the determinant of the Jacobian of F. We shall 
make use of the following result 

Lemma 2.2. Let F : Q, — >■ Q be a diffeomorphism. For u,v G H 1 (Q) m , let 
u = (F- l )*u := u o F G H 1 ^) 171 and v = (F _1 )*f G H 1 ^)" 1 . Then we 
have 

( 2 - 14 ) Q[A,B](u,v) = Q[F»A,F*B]{u,v). 
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Proof. It is verified directly by change of variables in integrations as follows 



The proof is completed. 



n±^A\^9xpdx a 1 } dy r 



a,P=l 



dv ,B 

{—u) v dy 

dy q J 



□ 



In the rest of this section, based on the above transformation properties, 
we give the construction of the two-dimensional cloaking devices, which we 
shall investigate in subsequent sections. We start our study by considering 
the cloaking of the unit central disc. In Section 4, we shall indicate how 
to extend our study to the general case. In the following, we denote by 
Bp the central disc of radius p and S p := dB p . Let {B2; A b , B b , /&} be the 
(regular) background/reference space, where f b is supported away from {0}. 
Here, we always assume that there is no eigenvalue problem for (|2.10p in 
the reference space, and hence there is a well-defined DtN operator on S2, 
namely Bb . . Consider the transformation F, defined by 



(2.15) F : 



B 2 \{0} -> B 2 \Bx 

x ^ (1 + ¥)R 

F blows up the origin in the reference space to Bi while maps B2\{0} 
to B 2 \Bi and keeps S2 fixed. We note that the blow-up transformation 
(|2.15|) has been extensively investigated for the design of three-dimensional 
cloaking devices in the literature (see jlj). Under the transformation F, 
the ambient reference medium in B2\{0} is then push-forwarded to form 
transformation medium in B 2 \Bi as follows 

(2.16) A c (y) = F*A b (x) and B c (y) = F*B b (x), yeB 2 \B 1 . 

In Section 3, we shall show that 

Theorem 2.3. {B 2 \Bi; A c , B c , f c } with f c = J'^F^yfi, and A c and B c 
given by A2.16\) is an invisibility cloaking device for the region Bi with respect 
to the reference space {B2; A b , B b , f b }. That is, for any extended object 



{B 2 ;A e ,B e ,f e } 



{B 2 \Bi;^ c ,S c ,/ c } in B 2 \Bi, 
{Bi;A a ,B a , f a } in B 1: 



where {Bi;^4 a ,i? a } is an arbitrary regular medium and f a G H 1 (Bi) m , 
have 



we 



( 2 - 17 ) A< A ev B e ,/ e - A A 6) B 6) /, 
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In the sequel, we conveniently define the push-forward of the source term 

as 

FJ = det(DF)(F~ 1 yf. 

So, the cloaking device in Theorem 12.31 is given by 

{B 3 \Bi; A c , B c , f c } = F,{B 2 \{0}; A b , B b , f b }. 

Next, we derive the explicit expressions of material parameters of the cloak- 
ing medium in Theorem 12.31 By (12.12p -( i2.13p and straightforward calcula- 
tions, we have for y E B 2 \Bi 

(2.18) A c {y) Jy\^U{y) ® I m A b (F^ (y))U(y) ® I m 

\y\ 

+ U(y) ® A b (F-\y))(I - II(y)) ® I m 



\y\ 



■(I-U(y))A b (F-\y))(I-U(y)) 



+ (I - U(y)) / m A(F- 1 (y))n(y) ® I m , 
(2.19) jBc (y)=l^lli ?fe (F- 1 (y)), 

where I m is the mxm identity matrix and H(y) : M 2 — > M 2 is the projection 
to the radial direction, defined by 

n(y)*=(VA^> 



2/1/ |y| 

i.e. n(y) is represented by the symmetric matrix |y|~ 2 yy r . For y G B2\Bi, 
we let y := y/|y| = [yi,y 2 ] T G Si and y x = [yfSy^F G s i be such that 
y ■ y 1 - = 0. Let l m G M m with each entry being 1, and set = y;l m G M m 
and = yi~l m G M m for Z = 1,2. Then, by straightforward calculations, 
together with the following facts 

U(y)y = y (I - U(y))y = 0, nfc/)^ = 0, (/ - 1%))^ = y\ 

we have 

E = E w^v))^. 

p,g=l P,9=l 



p,q=l P,9=l 



Since is a regular medium parameter, we have by (|2.4p that as 1 2/| ->■ 1+ 

2 2 

E \ AP oiy%)% ^ and E [^(f)^]*^ oo. 



P>9=1 P.9=l 
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Meanwhile, it is obvious that B c does not satisfy (12. 5p . That is, the algebraic 
conditions ()2.4|) and (|2.5p for a regular medium are violated by the cloak- 
ing medium {B2\Bij A c , -B c }, which exhibits both degeneracy and blow-up 
singularities as one approaches the cloaking interface Sf, i.e. from B2\Bi. 
We remark that this is in sharp difference from the study in [3] for three- 
dimensional cloaking devices where one would encounter only degeneracy 
singularities. 

3. Finite energy solutions for singular PDEs 

Consider the differential equation underlying the cloaking problem in The- 
orem E3J 

(3.1) V[A e ,B e ]U = fe on B 2 , u\ S2 = h. 

As we have seen in the end of Section 2, the cloaking medium {B2\Bi; A c , B c } 
is singular. One has to be careful in defining the meaning of a solution to 
the singular PDEs system (13. ip (see also Remark 13.91 for the necessity of 
introducing a suitable class of weak solutions to (|3.ip other than spatial 
^solutions). To that end, we define for 4>(y) € (f(B2) m , 



£[A,B\(4>) 



[ ( (A^d^yd^ + iB^yA dy 



1/2 



We note that if {B2; A, B} is a regular medium, then it is straightforward 
to verify 

( 3 - 2 ) £ [a,b](4>) ~ \\4>\\m(B 2 ) m - 

Here and in the sequel, for two relations Ri,R,2, Ri ~ R2 means that there 
exists two finite positive constants ci,C2 such that ciRi < R2 < C2R1. 
Also, in the following, for notational convenience, we shall frequently refer 
to Ri < R2 as Ri < C2R2- Next, one can verify directly that due to the 
blow-up singularity of A e on Sf, 



(3.3) 0e#(B 2 ) m and £ [AetBe] (0) < 00 iff 

00 Si 

Here and also in what follows, we make use of the standard polar coordinate 
5 J/2) | — > (t cos 9, r sin 9) in R 2 . Hence, we set 

(3.4) ^°°(B 2 r == G ^(B 2 ) m ; ^ g =0}, 
and 

(3.5) ^ °°(B2) m := ^°°(B 2 ) m n ^(B 2 ) m , 
which are closed subspaces of £ (B2)' m . Set 

W e {y) := 4(l ^ 7 l) I m for 1< \y\ < 2; I m for \y\ < 1. 

I y I 
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We note by using the fact B^ is a regular material parameter, together with 
(g39]) and (I23]l . that for <f> G ^°°(B 2 ) 

( 3 - 6 ) ^[4,W e ] (0) ~ % e ,Be] (</»)• 

On the other hand, it is easy to see that £[A e ,w e ](') defines a norm on 
^°°(B 2 ) m . Let 

(3.7) H{ A ^ Be] (B 2 r := cl{^ 0O (B 2 r;^ e] (.)} ) 

that is, the closure of the linear function space ^°°(B 2 ) m with respect to the 
singularly weighted Sobolev norm £[A e ,w e ](')- Clearly, one can consider the 
elements in Hi A B j(B 2 ) m as C m -valued measurable functions. Moreover, 
we have 

Lemma 3.1. The map 

cf>^D%J:= Yl A f^ 0G^°°(B 2 r, a = 1,2 

(5=1,2 

has a bounded extension 

(3.8) D% e : Hf Ae<Be] (B 2 r » M(B 2 ;C m ), 

where A4 (B 2 ; C m ) represents the space of complex C m -valued Borel measures 
on B 2 . Moreover, for u G Hy A B j(B 2 ) m , we have in the sense of Borel 
measures 

(3.9) {D%u)(S 1 ) = 0, a = 1,2. 

Proof. Let <f> G ^°°(B 2 ) m , g G C(B 2 ) m and 4> = F*4>,g = F*g G L°°(5 2 ) m . 
Then, it is straightforward to show that D A eft G L°°(B 2 ) m . Hence, it has 
that 

/ (D%cj ) ygdy= [ m^Tgdy 

= [ fz^fM mjr 9 d X 

Jbouo} \ tri, dx k d Xl y Bl 



On B 2 , one has dy a /dxk = O(-). This together with the facts that {B 2 ; Aj,, Bj,} 
and {Bi; A a ,B a } are regular, we use (|3.2p to further have 

mj)*9 dy\ 



B 



2 



< IHlH 1 (B 2 ) m ll5/ r ilL 2 (B 2 ) m + ||<^||iJ'l(B2) m ll5lU 2 (Bi)™ 

(3.10) < [f[A 6 A](^)+f[i4 a ,B ](^)]|bllc(B 2 )™dist(supp(^),Si; 

(3.11) < £ [AeijBe] (0)||5|| C (B 2 )mdist(supp(3),Si) 

(3.12) < 5[A e ,w e ](0)ll5llc(B 2 )'"dist(supp(5),Si) 



^ ll^llifi _ ,(B 2 )™llsllc(B 2 )™dist(supp(5),S 
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In the above inequalities, from (|3.10p to (|3.1ip . we have made use the fol- 
lowing facts by using Lemma 12.21 

(3.13) £[ Ab ,B b ](4>) = Q[A b ~B b ](4>,4>) = Q[F»A b -F»B h ]{4>,4>) = £[A e ,B R ](4>), 

whereas from (I3.11|) to (|3.12p . we have made use the equivalence (j3.6|) . This 
proves the bounded extension (13. 8p . Finally, ()3.9p follows by taking functions 
g supported in sufficiently small neighborhoods of Si. 

The proof is completed. □ 



Now, the solution to the singular system of PDEs (|3.ip is defined by the 
distributional duality as to find u G Hy A B j(B2) m such that u\s 2 = h G 

H 1 / 2 (S 2 ) m and 

(3.14) Q [Ae ,B e ](u,<f>) = </ e ,0>, V</> G ,%°°(B 2 ) m . 

Remark 3.2. Since the singularities of A e and B e are only attached to , 
we know that for any u G Hh B j(B2) m , u G i// oc (B2\Si) m . Therefore, for 

(|3.14p we have the well-defined u|s 2 = h G H 1 / 2 (S2) m and also a well-defined 
Dirichlet-to-Neumann map on S2 defined by 

(3.15) ^A M ( h ) = E ^Afdpu G ^- 1 / 2 (S 2 ) m , 

a,/3<l 

provided (|3. 14[) has a unique solution. 
Remark 3.3. As is known, 

f [ E {Afdpu)*d a u+{B e u)*u\ dy 

is the (generalized) energy of the system. Hence, the solution in (|3.14j) is 
physically meaningful in that it has finite energy. 

Next, we shall show the following on the solution of (|3.14p . 

Theorem 3.4. u G Hh B , (B 2 ) m is a solution to \3. 14\ ) if and only if 

v(x) = (i ? *(^| B2 \Bi)) e G H 1 (B 2 ) m with (i ? *(w|B 2 \Bi)) 6 denoting the exten- 
sion o/i ? *(u| B2 ^B 1 ) from B2\{0} to B2 (e.g., by setting it be 0), is a solution 
to 

(3.16) v [A b ,B b ]V = fb on B 2 , v\s 2 = h 
and w = u|bi G H l (Y±i) ni is a solution to 

(3-17) V[A a ,B a ] w = fa onBx, w\ Sl = c 0) 

where cq G C m is a constant vector determined by 

(3.18) I E "<»AZ%w dS = 0. 

1 a,/3<l 
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As a direct consequence of Theorem I3.4( we first give the proof of Theo- 
rem E31 

Proof of Theorem \2.3[ Let u G HK b j(B2) m be a solution to (|3.ip corre- 
sponding to the cloaking problem for {B2; A e , B e , f e }, whereas v G ff 1 (B2) m 
be a solution to (|3.16p corresponding to the differential equation in the ref- 
erence space {B2; Aj,, Bf,, ft,}. By Theorem 13,41 we have 

u Ib 2 \B! = (^ _1 )*(^Ib 2 \{o}), *Ib 2 \{0} = ^*Hb 2 \Bi)- 

Next, let n(t) G @(M) be a smooth real cut-off function such that < 
n{t) < 1 with n(t) = 1 for t < 4/3 and re(t) = for t > 5/3. For arbitrary 
4> G <f(B 2 ), set ^(y) := (1 - € ^°°(B 2 ). By Green's identity, 

we have 

Q[A ClJ B c ]KV0 = Q[ J 4e, J B e ](« ! V') 
= (/e,V'>B 2 + (A^ Ai/e (/ i ),7V'>S 2 

(3.19) =(/c,^)b 2 \B 1 + (^A e ,B e J e (h),^)s 2 . 

Then, by change of variables in integrations in (|3.19p . and using Lemma [2.21 
and, F|s 2 = Identity and f c = J^ 1 (F^ 1 )* fb, we further have 

(3.20) Q[A b ,B b ] (v, i>) = (fb,^B 2 + (A-A e ,B e ,f e ( h )^)s 2 , 
where = F*ip G $ (B2). By using Green's identity again, we know 

Q[A b ,B b ](v,-4>) = (fb,i>)B 2 + (^A b ,B b ,f b (h),7^)s 2 , 
which implies by (13.20|) that 

(Al )Be ,/ e (M,7^)s 2 = ( A A b ,B bJb (h),^)s 2 

and hence 

A A e ,B e J e ~ A A b ,B b ,f b - 

The proof is completed. 

□ 

We proceed to the proof of Theorem 13.41 We first derive two auxiliary 
lemmata characterizing the singularly weighted Sobolev space H} A B i(B 2 ). 

In the following, xn denotes the characteristic function for a set !!cl 2 . 
Lemma 3.5. Let c\ and C2 be two constant vectors in C m . Then 

(3.21) dXB, + c 2 Xb 2 \ Bi G H{ A ^ Be] (B 2 ) m . 

Proof. Since a constant function always belongs to B ^(B2) m , it suffices 
to show that ip(y) := XBi c £ H[A e B e ](-^ 2 ) m c G C m a constant vector. 
Let p G £F(R) be a cut-off function such that < p(t) < 1 with p(t) = 1 for 
t < 1/2 and p(t) = for i > 1. Then, define for e > 0, 

(3 - 22) ^ (y)= ( k y u % B B\s 

[ M ln(|y|-l) ) C ' yGB 2 \Bi. 
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Obviously, <p E G ^°°(B 2 ) m . Next, we shall show 
(3.23) lim \\tp e - ip\\ H} 



,(B 2 )* 



0. 



which then implies that G i?^ R 1 (B2) m . In fact, we have 

/ T3 „ \ 771 *- 

[A e ,S e 



,(B 2 ) m = £ [Ae,We}^ £ ~ 



(3.24) 



[A e ,B. 



B 2 \Bi 



By using the explicit expression for B c in (|2,19p together with the fact the 
B{> is a regular material parameter, it is straightforwardly shown that 

2 



(3.25) 



B 2 \B! 



(B c ip)*ip dy 



lne 



'B 2 \B 

< Voi(r 

where 

T £ := {y G M 2 ; 1 < |y| < 1 + e}. 
Next, by direct calculations, we have 



dy 



H\y\ - 1) 

-r as £-f0 + 



(3.26) V p( 



lne 



M\y\ - 1) 



lne 



lne 



M\y\-1)> |ln(|y|-l)P \y\-V 



where y :- 



(3.27) 



G Si. By using (LTTSI) and (EHf^ . we further have 



lne 



ln(|y|-l) 



)®c A* c Vp( 



lne 



lne 



1 2 



ln(|y|-l)' 



In el 



ln(|y| - 1) 
1 



|ln(|y|-i)| 4 |y|-i 



(A b (y®c))*(y®c)] 



By using (|3.27p . we can deduce 

/ Yl (Affy<Pe)*da<Pe dy 



(3.28) 



< 



^B 2 \Bi 

L 



|Ine| : 



lne 



MM-i) 



c A* c (Vp( 



lne 



ln(|y|-l)' 
as e^0+. 



c dy 



|ln(|?/|-l)| 4 M-l 

Combining (|5^]> . ([5I25]| and (f3728]) . we have (pH3j) . The proof is completed. 

□ 

Lemma 3.6. Let u be a measurable C m -valued function on B 2 . Then u G 
Hy A B j(B 2 ) m and only if the following two conditions hold: 
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(i) w := u|bi G fl^Bi)" 1 and 
(3.29) 7~w := u;| s - = constant, 

where 7" is the trace operator on ST, i.e. as one approaches S± 
from Bi. 
(ii) v := u|b 2 \Bi satisfies 



(3.30) 



(3.31) 



2 {Afd p v)*d a v dy 

B2 \ Bl ay9<l 



< OO, 



B 2 \B a 



(B c v)*v dy 



< 00 



and 



■ Vv. 



Pl S7 



0, p=l,2,...,m, 



where eg is the unit angular directional vector on the sphere S\. 



Proof. First, we show conditions (i) and (ii) are necessary for u 6 
Clearly, we have 

£[A a ,B a ](w)<\\u\\ H l AeBe]{B2)m , 

which together with the fact {Bi;A a ,B a } is regular implies that 



[A e ,B e 



fB 5 



lk||fli(Bi)» H U ll^ eiBe] (B 3 ) m ' 
Let {(/> n }n=i c ^°°(B 2 ) m be such that 



i.e., w e H 1 ^] 



u \\h} 



as n — > 00. 



[A e , Se ,( B ^ 

Then, we obviously have 

||<?Mbi - ^IIh^Bi)" 1 -> as w -> 0. 

Hence, 



|7 0n -7 w lltfi/2( Sl )* 



< 



|Bi — w \\H 1 (Bi) m ~~ ^ as n — )• 0. 



Noting dcfrn/dOls-L = 0, we know 7 n are constants independent of the 
angular variable 6. Therefore, j~w is constant. Whereas for v := u| B2 \g 1 , 
it is trivial to see that (I3.30P holds. To see (|3.3ip . we let v = F*v. Using 



the change of variables in (I3.30p . we have 

I / (A^dp)*d a v dx\ < 00 and | / (Bbv)*v dx\ < 00. 



a,B<l 



Since the reference space {B 2 ; A^, B^} is regular, we further have by using 
(ED and (1231) 



f EM ! *<» ««• / 

Jb 2 \{0} q<1 Jb 



\v\ 2 dx < 00. 



B 2 \{0} 
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Set % = VS p |b 3 \{0}, V = 1,2,..., m. We have % G L 2 (B 2 \{0}) 2xl . 
Extending ^ and v on {0} (e.g., by setting to be zero), and using the same 
notations for the extensions, we have 

* p G (L 2 (B 2 )) 2xl and v p G L 2 (B 2 ). 

The difference Vv p — ^ p belongs to i?~ 1 (B 2 ) and it is supported on {0}. 
Since no non-zero distribution in i/ _1 (B 2 ) is supported on {0}, we see Vv p = 
Hence, v p G H 1 (B 2 ) and therefore v G H l (W) m . 
Let the Fourier decomposition of v p , 1 < p < m, be given by 

(3.32) v p (r,9) = -}=J24 n) ( r y ne - 
In [2], it is proved 

(3.33) {)W(o) = 0, n^O. 

Formally, we write v p (0, 9) = v p °\o). Since v p x B , 2 \^l = {F^Yvp, we have 
V 1< r < 2, Up (r, 0) = -L ^ t> p n) (r)e me with v p n) (0 = v p n) (2(r - 1)). 
By (|3.33p . we have 



v 



P 



1) =«W(0) =0 for n t^O, 



which implies 

(3-34) -%M) 



r 5(9 



= 0. 

r=l 

Next, we show that (i) and (ii) are also sufficient conditions for a mea- 
surable function u on B 2 to belong to H^ A B j(B 2 ) m . We first assume that 
u Ib 2 \Bi = 0. By (i), let Co be a constant vector such that u\ s - = cq. Accord- 
ing to Lemma [331 it suffices to show that u := u — coXBi G Hh b j(Bi) m . 
Clearly, 

w := UXB! =w- CoXB! G ffd(Bi) m . 
Hence, there exists {ip n }^ =1 C £t(Bi) m such that 

W'Pn — wWffi^^m —7-0 as n — > oo. 

Let 4> n G ^°°(B 2 ) m be such that n | B i = </?n an d 0n|B 2 \Bi = 0- Using the 
fact that {Bi; A a , B a } is regular and u|b 2 \Bi = 0> we have 

" & % e>Be] (B 2 )™ =^ a ,B a ](^n " W) 

< ||v3 n - ^||^i {Bl) m -> as n -> oo, 

which implies that u G il?^ B ^(Bi) m . 

Now, let u be a measurable function satisfying (i) and (ii). As is shown 
above, u\b 1 G b ^ (B 2 ) m , one only needs to show that u — tiXBi G 



TWO-DIMENSIONAL INVISIBILITY CLOAKING 



15 



H[A e B e ] (-^2) m ; that is, u vanishes inside Bi. By (|3.3ip . we know u\ s + is 
constant independent of the angular variable 9. Using similar argument 
as earlier by substraction of a Heaviside function from u together with 
Lemma |3, 5\ we can further assume that u vanishes on Bi. Moreover, with- 
out loss of generality, we can also assume that u vanishes near S2. Let 
v := F * (w|b 2 \Bi) = F* v - According to our earlier argument, v G H (B2) m 
(note here we identify v and its extension on {0}). Since v vanishes near S2 
and {0} is a (2, l)-polar set, there are (f> n G ^(B2\{0}) m such that 

\\4>n ~ v\\m(B 2 ) m ~> as n -> 00. 
Let (j) n G J°°(B 2 ) m be such that 

0w|b 2 \Bi = (F^Y&n an d 0„| Bl = 0. 
Using u|g = and the fact that {B2; A^, Bj,} is regular, we have 

Un - uWHX^^m < f [B2 \B i; A c ,B c ](^n " «) 

^ l|0n - £||tfi(B 2 ) m as n 00, 

which implies that u G ^?u e £ e ](B2) m - 

The proof is completed. □ 

Now, we are ready to present the proof of Theorem 13.41 

Proof of Theorem \3.4\ We first show that if u G H^ A B j(B2) m is a solution 
to (13. lj) and (|3.14|) . then we must have the decoupled problems (|3.16p and 
pTTj) . 

Let v := u|b 2 \Bi an d ^ = In the following, as in the proof of 

Lemma l3.61 we identify v and its ^-extension from B2\{0} to B2. Clearly, 
v satisfies (|3.16p . On the other hand, let w := u|b 15 then is a constant 
vector independent of the angular variable by Lemma 13.61 So, it is trivial 
to see that w satisfies (13.17j) . Hence, we only need to show (13.18j) . whereas 
the determination of Co from (13.18P will be discussed in the subsequent 
Theorem 13.71 

Set 

£+ := {y G R 2 ; 1< |y| < 1 + e/2} , S" := {y G R 2 ; 1 - e/2 < |y| < 1}, 
and 

T+ := {y G R 2 ; \y\ = l + e/2} , T" := {y G R 2 ; \y\ = l- e/2}. 
We first show for 92 G ^ °°(B 2 ) m 



(3.35) lim 



J2j£ \a,/3<l 



dpu)*d a <p - u 2 (B e u)*ip - f*tp dy 
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In fact, by using the expression of (A e , B e ) in and the change of variables 
in integrations, together with f e being J^ 1 (F^ 1 )* fb in T,f and v being the 
i? 1 -extension of F*v from B2\{0} to B2, we have 

/ + ( E (Afd p uyd a y-u> 2 {B e uy<p-f> ) dy 
[ I {A?dpv)*d a <p - u; 2 (B b vr^ - / 6 V ) dx 



^ (J|£||ffi(B e )™ + ll/&llj?-i(B s:r ) ll^llfTHBe)™ ase^0 + , 

where <fi G if 1 (B2) m is the continuous extension of F* (0|b 2 \Bi) from B2\{0} 
to B2. On the other hand, noting u|b! = w G if 1 (Bi) m and / q |bi G 
if _1 (Bi) m , it is straightforward to verify 



(3.37) 



lim 

£-S-0+ 



: lim 

E-S-0+ 



^ E (^0%™)*^ - u 2 (B e w)^ - /> j dy 



0. 



Next, by (13. 14ft and using (|3.36j) and (13.3Tf) . together with integration by 
parts, we have 

(3.38) 







I I (Afdpuyd a <p-L0 2 (B e tpy<p-r e A dy 

2 \a,l3<l ) 



E (^V)*3«¥> " ^ 2 (^)V - / e > ] dy 

B2 \ Sl Va,/3<1 



e^0+ JB 2 \(S+USiUS e -) 



,a,/5<l 



: lim 



( J] l / 8 4V^(!')+ /_-( E VaAfdpuycpdS(y) 



o,/3<l 



a,/3<l 



where, by a bit abuse of notations, v denotes the exterior unit normal to 
respective domain, T,f and S~. Next, we estimate the integral in the right 
hand side of the last equation in (13, 38ft , In the sequel, we denote by v the 
exterior unit normal vector to the domain B e . Using the change of variables 
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in integration and the fact that v = -F*(w|b 2 \Bi) belongs to H (B2) m and 
satisfies (|3.16p . we have 



( J2 v*AfdpuyydS{y) 

a,B<l 



x; 



(3.39) 



dB, 



( v a Afdpv)*<pdS(x) 



a,B<l 



E \a,l3<l ) 

^ (ll^lli?i(B e ) m + II/6|Ih-i(b £ )») II^II^hb^)" 1 -> as e -> 0+ - 

By (EOHjl and (EOOjl . we see 

/ ( V Va AZ%w)*<pdS(y) = lim / ( V v a A<*/d p w)*ipdS{y) = 

which implies by the fact <p\si could be an arbitrary constant vector from 
C m that 

f ^Afdpw dS{y) = 0. 

Now, let v G if 1 (B2) m and w G i? 1 (Bi) m be solutions, respectively 
to (|3J6j) and (I3TT7D . Set u be (F" 1 )*^ on B 2 \Bi and be w on B b and 
extend it to B2 by setting it be zero on Si. By Lemma 13.61 we see that 
u G H^ A B j(B2) m . Moreover, it is readily seen that one also has 



(3.40) 

and 

(3.41) 



V[A C ,B C ]U = f c on B 2 \Bi, u\ S2 = h, 



v [A a ,B a \u = fa on Bi 



By Lemma 13. 1|. together with (|3.40p and (I3.41|) . we have for any ip G 

^ °°(B 2 r 

(3.42) 

Q[A e ,B e ](u,f) ~ (fe,<f)n 

= lhn f ( V (Afdpu)*d a <p - u\B e ^)*v - / e *0 I dy 

^0+ Jb 2 \(E+US 1 US £ -) I 



lim + I / + "( E ^Afdpu)*<pdS(y) + [ -( £ 



u a AZ%u)*<pdS(y) 
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Again, using the estimate in (|3.39p . we know 

(3.43) lim f ( V u a Afdsu)*ip dS(y) = 

e^0+ Jy+ 

E at,P<l 

On the other hand, by noting u = w on Bi, and using (|3.18p and <p\si is a 
constant vector in C m , we further have 

(3.44) lim f ( V v a Afdf}u)*<p dS(y) = 
Finally, by we have 

G[^,fl e ](«,y) = (fe,<Pfo, G ^o°°(B2). 

That is, u G Be] (B 2 ) m is a solution to (jBTTIj) . 

The proof is complete. □ 

In the rest of this section, we study the interior problem (|3.17p ~ (|3.18p . 
To that end, we introduce the following closed subspace of H (Bi) m , 

(3.45) W := {g G H x (&x) m ; 7 #| Sl = constant}. 
Then, ()3. 17j) ( 13. 18|) is weakly formulated as 

(3.46) weW and Q [Aa>Ba] (w,g) = (J a ,g) Bl , V<? G 
The corresponding homogeneous problem is 

(3.47) wdW and Q [Aa ,B a ](w, g) = 0, Vg G W, 
and by noting (12. 3p . its adjoint problem is 

(3.48) v £ W and Q[a,b*]O>»0) = °. V 5 G >T. 

Theorem 3.7. Let W denote the set of solutions to flff.^7| ). Then, either 
(i) W = {0}; or (ii) dim W = n for some finite n > 1. In £/ie case (i), the 
problem j3.^7[ ) is uniquely solvable. Whereas in case (ii), the homogeneous 
adjoint problem flff.^ffp a/so /ias exactly n linearly independent solutions, say 
vi,v%, . . . ,v n G W, and the inhomogeneous problem ft3.46\ ) is solvable iff 

(3.49) (v p , f a ) = /or 1 < p < n. 

Proof. Consider the operator C : W — >■ 3^* determined by Q[A a ,_B a ] i n the 
standard way as 

= G[^,,B ](fll,fl!i). 

Let = L 2 (Bi) m act as the pivot space. Clearly, the inclusion W C is 
compact. So, £ is Fredholm with index 0. Furthermore, each distribution 
f a G F _1 (Bi) m gives rise to a unique functional F a G defined by 

Fa (5) = {fa,g)'B 1 for 5 € Thus, the equation (|3.46p is equivalent to 

(3.50) Cw = F tt . 

The celebrated Fredholm theory applied to (|3.50p gives the desired results 
in the theorem. □ 
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Remark 3.8. In order to guarantee the existence of solutions to (13. ip . we 
impose the following compatibility condition on the interior source/sink f a , 

(3.51) («,/«> = 0, 

where v is any solution to (|3.48[) . 

Concerning the solution to the cloaking problem (I3.ip . we have the fol- 
lowing qualitative observations 

Remark 3.9. Let u 6 Hh B j(B2) m be a solution to (|3.ip . Set 7 + and 7" 
be the two sided trace operators on the cloaking Si. Then, by f|3.38[> . we see 



[ Y, v<*Afdpu] dS 



1 ot,P<l 

= [ [7 + ( £ "«^V) - 7 "( Y "«Afdf}U)] dS 
1 Oi,/3<l ct,P<l 

= 

On the other hand, we see that generically one has 

[u]\si '■= l + u — 7~n = constant vector 7^ 0, 

since otherwise we would have an over-determined system (I3.17p - (I3.18|) . 
This observation also encompasses the necessity of introducing the finite 
energy solutions other than spatial .ff ^solutions to the singular cloaking 
problem (|3.ip . 

4. General invisibility cloaking 

Up till now, we have considered invisibility cloaking where the cloaked 
region is the unit disc Bi, and the cloaking layer is the annulus B2\Bi. In 
this section, we extend our study to the general cloaking. To that end, let 
G : B2 — > £1 be an orientation-preserving and bi-Lipschitz mapping and let 
D = G(Bi). Then, for F in (1215]! . let 

K = GoFoG^ 1 : tt^Q, 

which blows up the point G({0}) to D within while keeps <9f2 fixed. Now, 
we shall show 

Theorem 4.1. Let {Q; At,, Bb, fb} with it supported away from G({0}) be 
the regular reference/background space. Then 

{n\D; A c , B c , f c } = K,{fi\G({0}); A b , B b , f b } 

is an invisibility cloaking device for the region D with respect to the reference 
space {17; Ab, Bb, fb}- That is, for any extended object 



{ft;A e ,B e ,f e } 



j{ftVD;A c ,B c ,f c } m Q\D, 
\{D;A a ,B a ,f a } m D, 
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where {D; A a ,B a } is an arbitrary regular medium and f a G H" 1 (D) m , we 
have 

Proof. Let 

(4.2) {B 2 ;A b ,B b ,f b } := (G^ 1 )*^; A b , B b , f b }, 

and 

{B 2 \Bi; A c , B C) / c } =(F o G^UnXGdO}); A b , B b , f b } 
=i ? 4B 2 \{0};A,B 6 ,/ 6 }. 

Then, 

A c ,B c ,f c } = G»{B 2 \Bi; ^Bc/c}- 
Since G is orientation-preserving and bi-Lipschitz, we know {B2; A b ,B b , f b } 
is a regular space and, A C ,B C have same singularities on the cloaking in- 
terface dD + as of A c , B c on S~[. In order to stick close to our discussion in 
Sections 2 and 3, we introduce 

{Bi; A a , B a , f a } = A a , B a , f a }, 

and define {B2; A e , B e , f e } correspondingly. Now we introduce the singulary 
weighted Sobolev space H}\ Be f ^(Q) m similarly to Hh B ^(B 2 ) m . The 
solution to the cloaking problem, 

(4.4) ^[A e ,B e ]U = f e on Q, u\ m = h 

can be defined in a similar manner as in (|3.14p . It is readily seen that 
u G B[A e ,B e ]( fi ) m is a solution to {44} iff u := G*u G Hf A ^ Be] (B 2 ) is a 
solution to 

(4.5) V[ AesBe ]U = f e on B 2 , u|s 2 = h 

where h = (G|s 2 )*h. By Theorem 13. 4^ we know v(x) = (F* (w|b 2 \Bi)) 6 e 
H l {B 2 ) m is satisfies 

(4.6) v \A b ,B b \v = h on B 2 , v\ S2 =h, 
whereas w = u\b 1 G /^(Bi)" 1 is a solution to 

(4.7) y [AaBa]W = f a onBi, w|si=e , 
with cq G <C m a constant vector determined by 



(4.8) 



I Y, v<*Afdpw dS = 0. 

•^ Sl a,B<l 



Therefore, we know v := (G" 1 )*^ = (K*(u| nv5 )) e G H 1 ^)™ satisfies 
(4-9) ^[A b ,B b ]V = f b on n, v\ m = h, 

whereas w =: (G~ 1 )*w = u\d G H l (D) m is a solution to 
^[A a ,Bal w = f a on A w Id = c o, 
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with Co G C m a constant vector determined by 

£ i/ a A.^w = 0. 

ay9<l 

Finally, by a similar argument to that for the proof of Theorem 12.31 one 
can show (|4.ip . The proof is completed. □ 

5. Discussion 

In this paper, we investigate the invisibility cloaking for a general system 
of PDEs in ]R 2 . The material parameters possess transformation properties, 
which allow the construction of cloaking media by singular push-forward 
transformations. The cloaking media are inevitably singular resulting in 
singular PDEs. Thus one needs to be careful in defining the meaning of 
solutions to the singular cloaking problems. In [3] the notion of finite energy 
solutions was introduced which is a meaningful physical condition because 
it relates to the energy of the waves. We follow their approach to treat more 
singular PDEs in the present paper. By studying the finite energy solution, 
we decouple the cloaking problem with one in the cloaking region having the 
same DtN operator as that in the background/reference space, and the other 
one in the cloaked region possessing some 'hidden' boundary condition on 
the interior cloaking interface. On the other hand, we would like to point out 
that this is not the unique way of defining solutions to the singular cloaking 
problems. An alternative treatment is provided in [7] by the authors to split 
the finite energy solutions space directly and hence decouple the cloaking 
problem into the cloaked region and cloaking region automatically. In order 
to achieve more physical insights into these solutions, one might relate them 
to solutions from regularized approximate cloaking. 
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